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Abstract 

We present in this letter a realistic construction of the coherent states for 
the Morse potential using the Klauder-Perelomov approach . We discuss the 
statistical properties of these states, by deducing the Q- and P-distribution 
functions. The thermal expectations for the quantum canonical ideal gas of 
the Morse oscillators are also calculated. 
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1 Introduction 



The one-dimensional Morse potential was first introduced as a useful model for 
diatomic molecules in 1929 [1]. Since its introduction, the Morse oscillator has 
proved very useful for various problems in diverse fields of physics and chemistry (di- 
atomic and polyatomic molecular systems, quantum chemistry, spectroscopy, chem- 
ical bonds) [2, 3, 4, 5, 6, 7] (and the references therein). 

On the other hand, since the pioneer papers of Glauber, Klauder and Perelomov 
about the coherent states of the one-dimensional harmonic oscillator [8], [9], and 
their different applications in physics [10], this approach has been paid attention 
also for other potentials. The three way to define the coherent states for the one- 
dimensional harmonic oscillator (i.e.: 1. as states which minimize the uncertainty 
relations; 2. as states which diagonalize the lowering operator and 3. as states 
obtained by action of the displacement operator on the bound state), which in this 
case are convergent (lead to the same results), for an another anharmonic potential 
lead to the three kind of states, generically named: Barut-Girardello, Klauder- 
Perelomov and Gazeau-Klauder [11], [12], [13]. 

Based on these papers, the coherent states was constructed for a series of poten- 
tials (see, [14], [15], [16], [17], [18]). 

For the pseudoharmonic oscillator the Barut-Girardello coherent states , as well 
as the photon-added Barut-Girardello coherent states was constructed [19], [20]. (see 
also, [21]). 

For the Morse oscillator, the construction of generalized coherent states was 
begin by the papers of [3]. Recently, Dong have constructed the coherent states 
based on the SU(2) realization for the Morse potential [22], Fakhri and Chenaghlou 
have constructed the Barut-Girardello coherent states for the Morse potential [23], 
while Roy and Roy , and, also, Popov have constructed and examined the properties 
of the Gazeau-Klauder coherent states for the same potential [24], [25]. So, the 
examination of the Klauder-Perelomov coherent states for the Morse potential, which 
is the main aim of this paper becomes a naturally task. 

This paper is organized as follows. In Sec. 2, after a brief recall of the main 
facts and properties of the Morse oscillator, we deduce the corresponding raising 
and lowering operators and we construct the Klauder-Perelomov coherent states for 
this oscillator. Section 3 is devoted to the examination of the statistical properties 
of these coherent states, by examining a quantum canonical ideal gas of Morse oscil- 
lators, in thermodynamical equilibrium with the reservoir, the concluding remarks 
will be given in Sec. 4. 

2 Klauder-Perelomov coherent states 

We shall begin by recalling some basic facts of the Morse potential, mainely the 
eigenstates and corresponding creation and annihilation operators. The hamiltonian 
for a one dimensional quantum particle evolving in the Morse potential is given by 
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where 



V\x) = (I + l) 2 - (21 + 3)e~ x + e' 2x (2) 
with I is an integer to simplify. The energies are given by 

E l n = (l + l) 2 -(l + l-n) 2 (3) 

where the quantum number n takes the finite sequence of values: 0, 1, I. 

The discrete spectrum of the Morse potential is finite and this fact is very im- 
portant in the construction of his coherent states as we will discuss in the sequel of 
this note. 

Because the ground state is zero energy it is well known that one can factorize 
the Hamiltonian as 

H l _ = X+Xr (4) 

where the operators Xf 1 (hermitian conjugated of each other, i.e (X^y = Xf) are 
given by 

X l ± = = F ^-(e--(i + l)). (5) 

From the supersymmetric quantum mechanics, it is well established that the 
supersymmetric partners 

K = X?X? = ~ + e~ 2 * - (2(1 + 1) ± l) e -* + (/ + I) 2 (6) 

have the same energy spectra, but different eigenstates. Indeed, it is easy to verify 
that 

H l + = H 1 - 1 + (21 + 1) (7) 

which implies that the eigenstates of H l + can be obtained from \ip l n ) corresponding 
to H l _ by substituting I by I — 1. The Schrodinger equations for Hi- are 



2 

where s — \. It is straightforward to check, for each I, the intervening relations 



H l ^n 2 ) = (l + lf-(l + -±s- n) 2 |^" 5±S ), (8) 



X?H l T = H l ± X?. (9) 

Due to the latter relations the operators Xf link the Hilbert spaces 7i l±s ~l 
and H lzfs ~i spanned by the bounded states of the Morse potential. They connect 
eigenfunctions with the same energy. Indeed one has 

Xl\€) = M eiam ~ n)+1) \^n-\) (10) 

KWn 1 ) = ^ie- iam ~ n)+1) \€ + i) (11) 



where a G R. The operators X l _ and X l + are not the ladder operators of the 
Hamiltonian H l _. To define the creation and annihilation operators for H l _, we 
define the operator 

^ = EI^)(C 1 I (12) 

rn 

satisfying U + U = 1 and iV'n-i) = ^Wn-i)- Using the transformation U, define a 
new pair of operators 

A l + = X l + U + , A l _ = UX l _ (13) 

Combining (10), (11) and (12), we get the actions of A l + and A l _ on the eigen- 

states (iV'n)} as 

^-l^> = V^ e<a(2(, " B)+1) l^n-l> (14) 

A l + m = V^e-^W^+Dl^) (15) 

The set of operators are the creation and annihilation operators of 

H l = A l + A l _ = X l + X l _. They satisfy the following commutation relations 

[A l _,A l + ] = -2N +{2l + l) (16) 

where the operator N is given by 

N\ti)=n\H). (17) 

The algebra generated by A l _ , A l + and N provides the appropriate tool to build 
up the coherent states for the Morse potential. 

The so-called Klauder-Perelomov coherent states for quantum system embedded 
in the Morse potential are defined 

\z,a) = exp(zA l + - zA l _)\4> l ) (18) 

Expanding the displacement operator and using the actions of the creation and 
annihilation operators, it comes 

k,a> = E^(N)e- fa ^K>. (19) 

n=0 

where 

00 (-\z\ 2 V 

4(W)=g^ijjTA'(„ + 1 ,i) (20) 
The quantities A which occur in the last formula are give by 

+ U) = (2l '\ nV E < E *4- E H (2i) 

(zi t- i u;. ii=1 i2=1 ._ 1 
with A'(n + 1, 0) = ^l+i-n)! • They satisfy the following reccurence equation 
A'(n + l,j) = yf2n(l + 1) - n?A l (n,j) + V2nl + 21 + 1- n 2 A l (n + 2, j - 1). (22) 
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Setting 



Jk(\z\) = \z\ n * 



(2J + l-n)! rI 



\ n\(2l + l)\ 
we get the first order differential equation 

dJ l n (\z\) 



4(1*1), (23) 



f^M) ~ (2ni + 2/ + 1 - n 2 )J l n+l {\z\). (24) 



d\z 

The solution of this equation takes the simple form 

1 



4(|*|) = ^(cosflzl))'-^ tg(\z\)r, (25) 
and the Morse coherent states rewrite as 

|Z,„> = (1 + |Z|T* E £^r,-+*M> (26) 

where Z = A They have the property of strong continuity in the label space 

and completeness in the sense that there exists a positive measure such that they 
solve the resolution to identity. The appropriate form of this resolution is 

/ dn(Z,Z)\Z,a)(Z,a\ = £ |<>«|. (27) 

J n=0 

To determine the measure dfi(Z, Z), we assume his isotropy and we set 

dfi{Z, Z) = (l + \Z\ 2 ) l h{\Z\ 2 )\Z\d\Z\d0/ 71 (28) 

with Z = \Z\e ld . Substituting Eq. (28) into Eq.(27), we obtain the following sum 
which should be satisfied by the function h(x = \Z\ 2 )) 

f°° t. w s , — n)\ , , 

/ x n h{x)dx = — ^ — ^. (29) 

The inverse Mellin transform gives 

This result can be obtained also by using the definition of Meijer's G-function 
and the Mellin inversion theorem [26]. 

Due to phase introduced in the definitions of annihilation and creation operators, 
the obtained coherent states are temporally stable. Indeed 

e-* tHi \Z,a) = \Z,a + t). (31) 

The physical utility of Klauder-Perelomov coherent sates for the Morse potential 
(26) in different applications consists in the calculation of the expectation (mean) 
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values of a certain physical observable A which characterizes the quantum system 
embedded in the Morse potential, with respect to \Z,a >: 



< Z, a\A\Z, a >=< A > z = (32) 
= (1 + l*| VjC^ (I) (^rV^-Ht) < MAM > (33) 

where we have used the binomial coefficients [27]: 

(')= " = T -^l . (34) 

[nj n\(l-n)\ F(n + l)T(l + 1 - n) 1 ; 

From a practical point of view, the most important operators are the diagonal 
operators in the eigenstates-basis \ip l n > corresponding to H l _, so that 

< ^ l m \A\ip l n >= a n 5 mn . (35) 

If A = N s , where the operator iV is given by Eq. (17) and s is an integer, we 
have: 

<A> Z = M 1 f (')ny = /i 1 „ (x4-) (l+x) 1 , (36) 

where we have used the simplified notation: x = \Z\ 2 . 
For s = 1 and 2 we obtain, successively: 

< N > z = l-^-, (37) 
l+x 

O t i/i \ ( 



These expectation values are useful in order to calculate the second-order corre- 
lation function for the morse oscillator [28]: 

2 < N 2 > z - < N > z I - 1 

{9)z = — WN^y — = ^ <L (39) 

Moreover, the Mandel Q-parameter is [29] 

Q z =< N > z \(g 2 ) z - ll = < 0. (40) 

These quantities provide information about the inherent statistical properties 
of the Klauder-Perelomov coherent states of the Morse oscillator \Z,a >. These 
properties depend on the analytical expressions of the functions (39) and (40) as 
depending on the variable x = \Z\ 2 . Because of the structure of these functions, 
{g 2 )z and Q z can be evaluated analytically. Generally speaking, the states \Z,a > 
exhibit sub-Poissonian statistics for those values of x = \Z\ 2 for which Q z < (or 
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(g 2 )z < l)(antibunching effect), Poisson statistics for values for which Qz = (or 
(fl ,2 )z = 1) and supra-Poissonian statistics for values of Z for which Q z > (or 
(9 2 )z > 1) (bunching effect). As we see from the above equations, the Klauder- 
Perelomov coherent states of the Morse oscillator \Z,a > are sub-Poissonian, for all 
values of the variable \Z\ 2 . 

By putting A = H l _ and using the eigenvalue equation (3), we obtain the action 
identity: 



<*->*=7T^£(iW = ( 41 ) 



= 2(1 + 1) < N > z - < N 2 > z = 1(21 +l)-^—-l(l- 1) (^-) = f( x ). 

1 + x \l + xj 

We see that this is a certain function of the label x = \Z\ 2 and this fact char- 
acterizes all the coherent states corresponding to systems with a finite dimensional 
energy spectrum, as it was mentioned earlier [24], [25]. 



3 Statistical properties 

We consider a quantum gas of the Morse oscillators in thermodynamic equilibrium 
with the reservoir (the thermostat) at temperature T, which obeys the quantum 
canonical distribution. Also, we consider that the individual Morse oscillators are 
in such states which are labelled by number state vectors \ip l n >- The corresponding 
normalized density operator for a fixed /-parameter (or, equivalently, for a fixed 
Morse oscillator) is then 

A = ^Ee-^|^><<l, (42) 

Z l n=0 

where (5 = (ksT)' 1 , ks is Boltzmann's constant and Z% - the normalization constant, 
i.e. the partition function for a certain fixed parameter /: 

Z X = E (43) 

n=0 

In the above exponential it appear the dimensional Morse eigenenergies: 

4 = = hum - -^—n 2 (44) 

n 2(1 + 1) n 2(1 + 1) y ' 

where uo is the angular frequency for the Morse oscillator with a fixed parameter /. 

The Q-function (the Husimi's function), i.e. the diagonal elements of the density 
operator in the representation of coherent states, is 

< Z,a\ Pl \Z,a >= ^j^ e -^( l \x n . (45) 
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It is not difficult to verify that the normalization condition of the density operator 
is accomplished: 



Trpi = J dfj,(Z,Z) < Z, a\pi\Z, a >= 1 
where we have used the following integral [27], [30]: 

,r(a)r(7 - a) 



dxx a ~\x + z)~~< = z a ^B(a, 7 - a) = z a ~ 



r(7) 



(46) 



(47) 



If we consider the dimensional Morse Hamiltonian, and the dimensional eigen- 
values of energy, instead of H l _ and E l n , then the corresponding energy exponential 
can be written as follows: 



(3e l n = fifiujn - (3- 



n = An — Bn 



(48) 



2(1 + 1) 

For most of the diatomic molecules B <C A. So, the limits of the parameter 
2(1 + 1) are very large, e.g. 37.1586 for H 2 molecule, i.e. for a "light" molecule and 
348.78 for I 2 , a "heavy" molecule [31]. As a consequence, the quantity B can be 
regarded as a perturbation constant and the energy exponential can be expanded in 
the power series as follows: 



-pe l „ 



OO K?fc 
-An \ " m 2k 



E 



n 



E 



B k 



k=0 nm k=0 

We can also use the following operator identity: 



d \ 
dAJ 



2k 



-A 



(49) 



E 

k=0 



dA 



2k 



exp 



B 



'A' 
dA. 



(50) 



Taking into account this ansatz, after the straightforward calculations, we can 
write the Q-function in the following manner: 



< Z, a\pi\Z, a >= — exp 
Zi 



B 



dA 



1 + xe A * 



1 



x 



(51) 



Because the Klauder-Perelomov coherent states for the Morse oscillator \Z,a> 
form an overcomplete set of states, they may be used as a basis set despite the fact 
that they are not orthogonal. Let us perform the diagonal expansion of the density 
operator pi in the coherent states basis: 



Pi 



jJdp(Z,Z)\Z,a>P l (\Z\ 2 )<Z, 



a\. 



(52) 



In order to find the quasi-probability distribution function Pi(\Z\ 2 ) (or the P- 
function) from the above diagonal expansion, we observe that the equation 



< f\Pi\9 >=jf dp(Z,Z) < f\Z,a> Pt(\Z\ 2 ) < Z,a\g > 



(53) 
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must be fulfilled for any arbitrary vectors < f\ and \g > from the Hilbert space (or, 
for any vectors from the basis \Z,a > or \ip l n >). 
The left-hand side of this equation is 

LHS = ^Y. e_/3d < f\€ >< €\9 >, (54) 

n=0 

while, after the angular integration 

± jyee-^-^ 9 = 5 nm , (55) 

the right-hand side becomes 



1 1 / ]\ fOO 1 

RHS =zM") <M>< *' 9>2(,+1) i ™ 2n+ ' jTTWY^ mZ?) - 

(56) 

After the variable change x = \Z\ 2 we must have: 



Jo 



°° dxx «- * » w = «-*-Lf'V = e -* 1 r(» + i)r(i + i-») 



/o (i + x)<+ 2 lv y z + iw z + i r(z + i) 

(57) 

where Pi(x) is an unknown function. In order to determine it, we extend the previous 
ansatz referring to the energy exponential [25]. 

It is obvious that the P-function also depends on the quantities A and B, besides 
the variable x. This leads to the idea that Pi(x) = Pi(x, A, B) can also be expanded 
in a power series similarly to the energy exponential (49) in the following manner: 



oo r>k 



B k 



k=0 



B=0 



00 B k ( d \ 2k 



where the function Xi(x, A) is to be determined. 

By inserting Eqs. (49) and (58) into the equation (57), performing the function 
change obtain: 



X t (x,A) 



1 



[l + xY+X&A) 



(59) 



r(/ + 2) 

as well as by extending the natural values of n to complex s such as n — s — 1, we 
get to the following Stieltjes moment problem: 



roc 1 

/ dxx^hAx.A) = e A — —T(s)T(l + 2-s). 
Jo (e- 4 )* 



(60) 



The solution of such a problem is [26] 
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h l (x,A) = e A T(l + 2) 



1 



1 + e A x) l + 2 ' 



Finally, the P-function is 



Pi(x) = exp 



B 



dA 



,A 



1 + x 
1 + e A x 



1+2 



(61) 



(62) 



where we have use the operator identity (50). 

The scalar product of two Klauder-Perelomov coherent states for the Morse os- 
cillator is 



< Z,a\Z',a >- 



(l + ZZ 



:i + \z\ 2 )Ui + \z'\ 2 ) l 2 



(63) 



From the trace condition (Trpi = 1) applied to the diagonal expansion (52) 
and the scalar product, it is not difficult to prove that the P-function satisfies the 
normalization condition: 



1 

~Z, 



dp(Z : Z)P l (\Z\ 2 ) = 1. 



(64) 



Let us verify this equation using the our obtained expression for the P-function 
(62). Thus, we have 



T 1 n n , f 2 - dOl r°° J 1 

1 = — (/ + 1) / / dx- — exp 

Z t y 'Jo it 2 Jo (1 + x) 2 1 



'a 



-A 



1 + x 
1 + e A x 



1+2 



= (65) 



hi + 1) exp 



e A dxil+xYil + expAx)' 
Jo 



The last integral is of the following kind [27], 



°° {ax + bf- 1 (adf - {bcf 



o (cx + dY+ l (3(ad-bc)(cdy 
This lead to the results: 



[cd>0; ad^bc- Re/3 > 1]. 

(66) 



1 = 4- ex p 



A\l+1 



1 ~ (e- 4 ) 
(1 - e A )(e A ) 1 _ 



— exp 



5> 



-A\n 



ln=0 



1. 



(67) 



So, our obtained expression for the P-function (62) is correct. 
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In this way, the diagonal representation of the normalized density operator of 
the Morse oscillator in the Klauder-Perelomov coherent states representation is 



Pi 



Y j MZ, Z) 



exp 



B 



'A" 
dA, 



-.A 



1 + 


Z 


2 


l + e A \ 


z 


2 



\Z,a><Z,a\. (68) 



Then the thermal expectation value (the thermal average) of an observable A 
concerning the Morse oscillator is given by 

<A> l =Tr{p l A) = jJ dii{Z^P x {\Z\ 2 ) < Z,a\A\Z,a > . (69) 

If the operator A is an integer power s of the number operator N, then, using 
Eqs. (69), (36, (47), as well as (49), we obtain the expected result: 



n e 



J l n=0 



(70) 



With this result we can define and calculate the thermal second-order correlation 
function (fi 1 ^) and the thermal Mandel parameter Q h i.e. the thermal analogue of 
the corresponding functions for the Klauder-Perelomov coherent states \Z, a > (see, 
Eqs. (39) and (40)): 



< N 2 >i — < N >i 
(<AT> ; ) 2 



1 + 



+ 



2 ' 



Q X =< N>1 [(^)) i] = _i- 



f^ 2 lnZ ; 



\nZ, 



(71) 



(72) 



dA ^ ^ l 

The difference between these results and the similar results for the Gazeau- 
Klauder quasi-coherent states of the Morse oscillator [25] (i.e. the absence of the 
term 1/2) is due definition of the Morse energy eigenstates (3), comparing by E n = 
hu(n + 1/2) - 2{l+1) \: +l/2) , (see, Eq. (3) in Ref. [25] ). 



4 Discussion and outlook 

In conclusion, we have given a construction of Klauder-Perelomov coherent states 
associated with the bound states of the Morse potential. These states possess all prop- 
erties of the usual coherent states as continuity, over completion, non-orthogonality 
and temporal stability, being defined on the entire complex space of the variable Z. 

In order to examine the statistical properties of these states, we have concen- 
trated our attention on a quantum system which consists on the quantum canoni- 
cal ideal gas of Morse oscillators, in thermodynamic equilibrium with the reservoir 
(thermostat). Firstly, we have obtained the expressions of the diagonal elements 
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of the density operator in the Klauder-Perelomov coherent states (the Husimi's Q- 
function), as well as the diagonal representation of the density operator and the 
corresponding P-function. With these functions we are able to express the ther- 
mal expectations (thermal averages) for the physical quantities which characterize 
the Morse oscillator quantum system. The second-order correlation function and 
the Mandel parameter, as well as their corresponding thermal pairs, offer the in- 
formation on the character of the coherent states, versus the Poisson distribution 
function. 

Besides the implicit construction of the Klauder-Perelomov coherent states for 
the Morse potential and examination of their properties, the main results of this 
paper are the derivation of the integration measure and the deduction of the ex- 
pressions for Q- and P-functions and, consequently, the thermal averages calculated 
with these functions. In our opinion, the above obtained results seems to be entire 
new, because, to our knowledge, these results have not yet been published in the 
literature. 

By using the general expression for thermal averages one can easy be deduced 
some interesting concrete thermal averages for the Morse oscillators quantum canon- 
ical ideal gas as: free and internal energy, entropy and molar heat capacity at the 
constant volume and so on. 

Acknowledgements: M.D is grateful to Max Planck Institute for the Physics 
of Complex systems for its support . 
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